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L AMENDMENTS TO THE CLAIMS : 

Kindly amend claims 1, 4, 7, 9, 11, 14 and 15 as follows. 

The following Listing of Claims replaces all prior listing, or versions, of claims in the 
above-captioned application. 
Listing of Claims: 

1. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, directly using V-CAD data, comprising the steps of: 

(A) dividing external data into a plurality of cells having boundaries orthogonal to 
each other, the external data including boundary data of an object that contacts 
incompressible viscous fluid; 

(B) classifying the divided cells into an internal cell positioned inside or outside the 
object and a boundary cell including the boundary data; 

(C) determining cut points in ridges of the boundary cell on the basis of the boundary 

data; 

(D) determining a polygon connecting the cut points to be cell internal data for the 
boundary face; and 

(E) applying a cut cell finite volume method combined with a VOF method to a 
boundary of a flow field to analyze the flow field, wherein step (E) comprises the steps of 

i. applying a two-dimensional QUICK interpolation scheme to a convection 
term for space integral; 

ii. applying a central difference having precision of a degree of a second 
order to a diffusion term; 

iii. combining the convection term and the diffusion term, and applying an 
Adams-Bashforth method having precision of a degree of a second order to the 



-2- 



Patent Application Serial No. 10/540,861 
Attorney Docket No. ASAIN0165 

combined convection term and diffusion term for time marching; and 

iv. applying a Euler implicit method having precision of a degree of a first 

order to a pressure gradient term for time marching, 

wherein for a two-dimensional boundary cell, a governing equation in the finite 

volume method is expressed by a governing equation (7), 

jj^dV = - jjdiv(u ® u)dV - \\div(pI)dV + J- \\div{grad(u))dV ( 7 ) 

wherein corresponds to the convection convention term, 

- §div{pf)dV 

u corresponds to the pressure gradient term, and 

corresponds to the diffusion term , wherein £ designates velocity of 
flow of viscous fluid, V designates differential volume of the viscous fluid, pi designates 

pressure p of the viscous fluid along the I vector, and Re corresponds to a non-dimensional 
Reynolds number; and 

(F) outputting to an output device a result of the method for numerical analysis of the 
flow field of incompressible viscous fluid, wherein the output device prints, or displays, or 
prints and displays, the result . 
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3. (Cancelled) 
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4. (Currently Amended) A method for numerical analysis of a flow field of 

incompressible viscous fluid, according to claim 1, wherein the convection term, the pressure 

gradient term and the diffusion term in the governing equation of the finite volume method 

are expressed by the equations (8), (9) and (10), respectively, 



convection term: 

ljdiv(u <S> u)dV = <g(w ® u) HdS = £ ( 5 ® " "* y -. 

+AX(4 jW ^ /2(/+1/2 V,, ^ 4j-l U i-V2J-U2 V IJ-\ )]' ( 8 > 

+Ax(4 ^.v^v. 7 . - 4 v^v, J _ 1 )]j\only no - slip on wall 
pressure gradient term: 

fldiv{pi)dr = c$(piynds= X pJ- nss » 

v i.j 

= AAB ijPi+V2J -B,._ ]jPi _ ll2J -p p {B tJ -B„Ji$ ( 9 ) 

+Ax[A i j p i J+u2 - A i M p UJ _ U2 - p p {A uj - A LJ _,)y 

diffusion term: 

$div{grad(u))dV = <$grad(ii) ■ ndS = £ grad(u) m • ii5S m 

v ..} s i-s m=l " 5 

= [Ay(5, jgrad(u)l U2J tt zradW^ -(£,, -5 f _, .^^(m)*) 
+Ax(4.,.grarf(«)^ (10) 

+Ax(4, ; .^(v)^, 1/2 -4..,gmrf(v)>. 1/2 -(4 >; -4^)«wrf(v)J)]7 

, wherein 

Ax is grid width in the x direction and Ay is grid width in the y direction, and vectors i and 

j correspond to unit vectors in the x and y directions, n corresponds to a unit vector normal 

to the differential surface S, and wherein A^ t Kj - u Bn, and B/ _^ are fractions of area of grid, 
and wherein u designates velocity of flow in the x direction and v designates velocity of flow 



-4- 



Patent Application Serial No. 10/540,861 
Attorney Docket No. ASAIN0165 

in the y direction . 

5. (Previously Presented) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 1, wherein when a no-slip boundary 
condition is used for a solid boundary, a first integral is performed on the solid boundary with 
the convection term being zero, a value of a middle point P of a cut line segment is used as an 
average for the pressure gradient term and the diffusion term, and a space integral is 
performed with areas fractions being applied to all of the terms. 

6. (Previously Presented) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 1, wherein the boundary cell having the 
parameter smaller than a threshold value of VOF=0.01 is regarded as a complete solid, 

for the boundary cell having the parameter larger than the threshold value, a 
definition point for the parameter calculated in a cut cell is set at a center of the boundary 
cell, 

and a definition point for a parameter in a ridge is set at a center of a cell ridge, and a 
parameter at a middle point of a line segment is calculated by a linear interpolation. 

7. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 1, wherein a drag force in a flow direction 
and a lift force in a direction vertical to the flow, acting on the object, are expressed by 
equations (12) and (13) that respectively express drag force and lift force as follows, 
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drag force : 




only Cartesian 



(12) 



lift force: 



v v °y s s 



(13) 



only Cartesian 




, wherein Fx 

designates drag force in the x direction and F Y designates lift force in the y direction, and a 
designates a stress tensor at the surface of the object. . 

8. (Previously Presented) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 1 , wherein in fluid-structure interaction 
analysis accompanying a moving boundary, a fluid system and a structure system are 
separately analyzed for each predetermined time interval, and boundary conditions for the 
fluid system and the structure system are explicitly used. 

9. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 8, wherein in analysis on a forcibly vibrated 
circular cylinder, the circular cylinder is set as a one-mass-point and a one-degree-of- freedom 
system such that the circular cylinder is a solid structure elastically supported and vibrating in 
a direction vertical to the flow, 
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and v-directio n Y direction displacement of a center of the circular cylinder is given 
by equation (17), and a velocity boundary condition in the v-directio n Y direction for a 
surface of the circular cylinder is given by equation (18), 
y=Asin(2 tz f c t) ...(17) 
v w =A2 7i f c cos(2 n f c t) ...(18) 4 

wherein A designates amplitude of v-direction displacement v, k designates 
frequency of forced vibration, and t designates time . 

10. (Previously Presented) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 9, wherein movement velocity of the 
vibrating circular cylinder obtained by the equation (18) is changed to be given for each 
calculation time step for the velocity boundary condition on the flow field. 

11. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 8, wherein in analysis on self-induced 
vibration due to a vortex shedding from the circular cylinder, a vibration equation having a 
dimension is expressed by equation (19) or equation (20), using a one-mass-point and a one- 
degree-of- freedom dumper/spring model, 

m 7^ +c ^ +k y=\p u « DC L < 19 > 

^ H 4,kfJ^H2^y-fc L (20) 
dt dt Sc 

wherein m designates mass of the circular cylinder, c designates a viscous damping 
coefficient, k designates an elastic recovery coefficient of the spring model, p designates 
fluid density. Up designates velocity of the uniform flow, D designates a diameter of the 
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circular cylinder, and C L designates a lift force coefficient, 

wherein fo designates a characteristic frequency, h designates a non-dimensional 
damping coefficient, Sc designates a Scruton number, and y = Asin(2 iz f ct), 

wherein A designates amplitude of y-direction displacement y, f £ designates frequency 
of forced vibration, and t designates time . 

12. (Previously Presented) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 11, wherein the movement velocity of the 
vibrating circular cylinder calculated by equation (20) is changed to be given for each 
calculation time step for the velocity boundary condition on the flow field. 

13. (Original) A method for numerical analysis of a flow field of incompressible 
viscous fluid, according to claim 11, wherein initial displacement and initial velocity of the 
circular cylinder are set to be zero, the lift force is explicitly given by using a current value, 
and the vibration equation is integral by the Newmark's j3 method to obtain vibration 
displacement and vibration velocity of the circular cylinder. 

14. (Currently Amended) A device for numerical analysis of a flow field of 
incompressible viscous fluid, directly using V-CAD data, the device comprising: 

an input device for inputting external data including boundary data of an object that 
contacts incompressible viscous fluid; 

an external storage device for storing substantial data of shape data and physical 
property data integrated into each other, and a storage operational program for the substantial 
data; 

an internal storage device and a central processing device for executing the storage 
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operational program; and 

an output device for outputting a result of execution of the storage operational 
program; 

wherein the device for numerical analysis 

(A) divides the external data into a plurality of cells having boundaries 
orthogonal to each other; 

(B) classifies the divided cells into an internal cell positioned inside or outside 
the object and a boundary cell including the boundary data; 

(C) determines cut points in ridges of the boundary cell on the basis of the 
boundary data; 

(D) determines a polygon connecting the cut points to be cell internal data for 
the boundary face; and 

(E) applies a cut cell finite volume method combined with a VOF method to a 
boundary of a flow field to analyze the flow field, wherein when the device for 
numerical analysis applies the cut cell finite volume method combined with a VOF 
method to the boundary of the flow field to analyze the flow field, the device for 
numerical analysis operates to 

i. apply a two-dimensional QUICK interpolation scheme to a 
convection term for space integral; 

ii. apply a central difference having precision of a degree of a second 
order to a diffusion term; 

iii. combines the convection term and the diffusion term, and applies 
an Adams-Bashforth method having precision of a degree of a second order to 
the combined convection term and diffusion term for time marching; and 

iv. applies a Euler implicit method having precision of a degree of a 
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first order to a pressure gradient term for time marching, 

wherein for a two-dimensional boundary cell, a governing equation in the finite 

volume method is expressed by a governing equation (7), 

j&JV = - \\div(u <S> u W - \\div{p!)dV + ^- §div(grad(u)~)dV ( 7 ) 

wherein corresponds to the convection convention term, 



corresponds to the pressure gradient term, and 

corresponds to the diffusion term , wherein u designates velocity of 
flow of viscous fluid, V designates differential volume of the viscous fluid, pi designates 

pressure p of the viscous fluid along the I_ vector, and Re corresponds to a non-dimensional 
Reynolds number . 

15. (Currently Amended) A computer readable medium comprising a program 
stored thereon for numerical analysis of a flow field of incompressible viscous fluid, directly 
using V-CAD data, wherein the program causes a computer to perform the steps of: 

(A) dividing external data into a plurality of cells having boundaries orthogonal to 
each other, the external data including boundary data of an object that contacts 
incompressible viscous fluid; 

(B) classifying the divided cells into an internal cell positioned inside or outside the 
object and a boundary cell including the boundary data; 

(C) determining cut points in ridges of the boundary cell on the basis of the boundary 

data; 
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(D) determining a polygon connecting the cut points to be cell internal data for the 

boundary face; and 

applying a cut cell finite volume method combined with a VOF method to a boundary 
of a flow field to analyze the flow field, wherein step (E) comprises the steps of 

i. applying a two-dimensional QUICK interpolation scheme to a convection 
term for space integral; 

ii. applying a central difference having precision of a degree of a second 
order to a diffusion term; 

iii. combining the convection term and the diffusion term, and applying an 
Adams-Bashforth method having precision of a degree of a second order to the 
combined convection term and diffusion term for time marching; and 

iv. applying a Euler implicit method having precision of a degree of a first 
order to a pressure gradient term for time marching, 
wherein for a two-dimensional boundary cell, a governing equation in the finite 
volume method is expressed by a governing equation (7), 



jj^-dv = -jjdiv(u ®u)dV - \\div{pl)dV^^- \\div{grad{u)W 



(7) 



-§div{u<&u)dY 

wherein *' <' corresponds to the convectio n convention term, 

' corresponds to the pressure gradient term, and 

corresponds to the diffusion term , wherein ii designates velocity of 



Re 



flow of viscous fluid, V designates differential volume of the viscous fluid, pi designates 
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pressure p of the viscous fluid along the I vector, and Re corresponds to a non-dimensional 

Reynolds number; and 

(F) outputting to an output device a result of the method for numerical analysis of the 
flow field of incompressible viscous fluid, wherein the output device prints, or displays, or 
prints and displays, the result . 
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